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Recovering modified Newtonian dynamics by changing inertia 
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ABSTRACT 

Milgrom's modified Newtonian dynamics (MOND) has done a great job on accounting for the 
rotation curves of a variety of galaxies by assuming that Newtonian dynamics breaks down 
for low acceleration typically found in the galactic contexts. This breakdown of Newtonian 
dynamics may be a result of modified gravity or a manifest of modified inertia. The MOND 
phenomena are derived here based on three general assumptions: 1) Gravitational mass is 
conserved; 2) Inverse-square law is applicable at large distance; 3) Inertial mass depends on 
external gravitational fields. 

These assumptions not only recover the deep-MOND behaviour, the accelerating expan- 
sion of the universe is also a result of these assumptions. Then Lagrangian formulae are de- 
veloped and it is found that the assumed universal acceleration constant cir, is actually slowly 
varying by a factor no more than 4. This varying 'constant' is just enough to account for the 
mass-discrepancy presented in bright clusters. 
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1 INTRODUCTION 



The modified Newtonian dynamics (MOND), originally proposed by Milgrom ( 1983) as an alternative to the cold dark matter paradigm to 
acco unt for the rotation curves o f spiral galaxies, has extended its success to dwarfs, low surface brightness galaxies (LSB) and ellipticals 
(see ISanders & McGaugh 1 2002, for a review). When confronting with clusters, especially rich clusters, MOND shows some drawback. 
On the cluster scale, MOND still needs dark matter, which is what MOND was particularly devised to eliminate. To o vercome this dif- 
ficulty, neutrinos were speculated to be responsible (Sanders 2003: lAngus et al. 1 120071 : iGentile. Zhao & Famaev 1 12008). Neutrinos with 
mass ~ 2eV, marginally allowed by current most accurate neutrino mass measurement, contributing negligibly to galaxies' mass budget, 
could be dynamically significant in clusters of galaxies. Th ough this hypothesis is successful in some aspects, it is still controversial (e.g. 
Pointecouteau & Silkl2005l : lAngus, Famaev & Buote 1120081) . 

Despite this drawback, MOND has drawn much attention because of its impressive success compared with the standard cold dark 
matter paradigm, which is facing with some difficulties, especially on the galactic scale. MON D is a phenomenological theor y that may 
be interpreted in different ways. First of all, it may indicate a breakdown of Newtonian gravity dBekenstein & Milgrom I I 1984 ^ where the 
standard Poisson equation is replaced by V ■ [fj, (| V</j| /ao) V<^] = AitGp, and ao (~ 1.2 x 10 -8 cm s -2 ), introduce d by MOND, is 
a new acceleration constant below which dynamics and/or gravity become significantly non-Newtonian. I Milgrom" 1 2002 ) reviewed this 
interpretation which relates the gravity to a potential flow. Although this is a field need more investigation, the lack of profound physical 
foundation ma kes this inte rpretation less attractive. A second interpretation is that gravitational constant increases when accelerations are 
lower than ao dBekenstein 2004). This relativi stic extension of MOND can mimic MOND's behaviour at low acceleration extreme, but it is 
still a subject of debate (e.g. Reves et al. 1l2010h. 

A third interp retation is what | Milgrom 1 1983 ) proposed that the Newtonian dynamics may break down at low accelerations. Instead of 
the usual F = ma, Milgrom 1 19831) suggested a modified dynamics 



F = m g n (a/ao) a, 



(1) 



M(a;> 1) 



(j, (x <C 1) w x, 
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where m g is the gravitational mass of the body moving in the field. This relation is equivalent to 

I \l/2 

a » (a N ao) 



(2) 



in the deep-MOND regime, where is the acceleration derived from Newtonian dynamics. It is just this simple relation that works remark- 
ably well in reproducing the dynamics of a variety of galaxies with quite different morphologies and luminosities. 

The modified dynamics can be interpreted as a modification of inertia, which, if applied properly, may solve a lot of puzzles faced with 
modern physics. Modification o f inertia is no t a new idea since the time of Mach who ch allenged Newton's idea about inertia. Since the 
theorization of Unruh radiation i Unruh lH976h and the discovery of Casimir effect i Casimirl 194S ), their relation to inertia and gravity has 
frequently been speculated by several authors (e.g. Puthoff |[l989l : Haisch. Rueda & Puthoff II 1994 ). Though Casimir effect is a reality both 
theoretically and experimentally, what does it mean for and how to apply it to cosmology remains a subject of debate. Puthoff 's gravity, 
based on Casimir effect, has the flaw of missing experimental support. As for Unruh radiation, it is not clear what this radiation means for 
cosmology and whether it is related to inertia. 

In this paper, we propose a new approach, based on some general speculations, to modification of inertia. By this modification of inertia, 
Q is successfully reproduced. Lagrangian formulae are developed and its indications are discussed then. 



2 ASSUMPTIONS AND RESULTS 
2.1 Assumptions 

In this section, three assumptions are proposed: 1) Gravitational mass is conserved; 2) Inverse- square law is applicable at large distance; 3) 
Inertial mass depends on external fields. 

If we make a comparison between gravity and electromagnetic interaction, we immediately realise that gravitational mass is analogous 
to electric charge. It is the electric charge, as a source, who produces electric field. Electric charge cannot be created and destroyed, it is just 
a being. Yang-Mills gauge theory, the foundation for standard model in particle physics, is based on the conservation of charge. As a source 
of gravitational field, it is unphysical that gravitational mass is not conserved. 

In electromagnetic field, the inverse-square law is directly rel ated to the zero-mass-ness of photons . Coulomb's law has been tested from 
~ 2 x 10 10 m down to 10 -18 m, a magnitude span of 28 orders i Adelberger. Heckel & Nelson 2003). The astronomical tests of Newton's 



gravity has been mainly confined within solar system. The most accurate astronomical tests are lunar-laser-ranging studies of the lunar orbit, 
which do n ot show a deviation of gra vity from Newton's law. On galactic scale, a Yukawa-like gravity has been proposed by many a uthors 
(see, e.g. Sanders & McGaugh 2002, and references therein) to accou nt for the mass discre pancy in galaxies. However, Milgrom 1 1983h 



pointed out that this is inconsistent with the empirical Tully-Fisher law JTullv & Fisher II 19770 . This line of arguments indicate that gravity is 
not Yukawa-like from submillimeter scale to at least galactic scale. Large scale structure of the Universe favors a Newton's gravity law even 
on cosmic scale. Therefore it is quite safe to assume that inverse-square law is accurate on the scale of our interest in this paper. 

The equiv alence of gra vitational mass to inertial mass, upon which general relativity is based, has been tested experimentally with very 
high accuracy Iwilll l2009). But these tests are confined within solar system, no direct test is available on the galactic scale. Unlike the 



gravitational mass, inertial mass does not associate with any physical field. Physically, because of the association with gravitational field, 
gravitational mass should be conserved, but this is not necessarily the case for inertial mass. Inertia is one object's ability to keep its original 
state of motion. However, how does one object know its original motion state if no reference is available? 

Mach speculated that inertia is the result of the object's motion relative to the mean mass distribution of the Universe as a whole. In 
other words, inertia is meaningless if no mass is there other than the object itself. Let us consider another situation: If the mass is distributed 
uniformly throughout the Universe but one distinct object. Whatever the state of motion of this distinct object, the state of the Universe will 
keep the same. This is indicative of the dependence of inertia on the mass distribution of the Universe. 

It has been long that Milgrom found that external field plays a role on the internal dynamics of open clusters where the internal field 
is well below the critical acceleration ao and therefore should be in deep-MOND regime. However, the dynamics of these systems does 
not show any evidence of dark matter. Milgrom realised that external field, if significantly above the transitional acceleration between 
Newtonian dynamics and MOND, could play a role and make the dynamics of these open clusters Newtonian. In addition, according to 
Mach, gravitational field, being the result of the mass distribution of the Universe, therefore endows the object with inertia. This is slightly 
different from the original Mach's Principle. Other than simultaneously depending on the mass distribution of the whole Universe, here we 
speculate that inertial mass depends on the local gravitational field. This speculation is based both on Milgrom's realisation and on the present 
day basic belief that there cannot exist action-on-distance. 

The above lines of argument indicate that if no external gravitational field exists, the object could be in random states of motion, i.e. the 
inertial mass is zero. With the increase of external field, inertial mass increases accordingly. The climbing-up of inertia, however, does not 
continue infinitely. Having the external field, i.e. the gradient of the potential, as a reference, the object has a sense of its past motion state. 
The stronger the external field, the more sense it has about its original motion state. But if the external field is strong enough, increasing the 
field's strength will not increase its sense of past motion state because it just has enough "information" about its original state of motion. As 
a result, we set the inertial mass mi in the strong field limit to be its gravitational mass m g , as we know from the dynamics of solar system. 

However, in what a way the external gravitational field tells the moving particle of its past motion state? The only way that this can be 
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achieved is by the interaction of the particle's historical trajectory with the particle itself. That is to say, some present-day unrecognised field 
(referred to as S field hereafter) is released when the particle is in motion. If the S field is freely moving without any interruption, it will never 
interact with the particle who released it. In this case, the particle does not know its past motion state and is therefore inertialess. But if the S 
field it released is interrupted, the particle will acquire inertia because of its interaction with its prior history. 

This argument seems still to be somewhat speculative. Let us now discuss a real situation where an inertialess object could naturally 
acquire inertia as a result of its interaction with its prior history. It is well known that crystal is never perfect and full of varying kinds of 
defects, e.g. line defects and point defects. Line defects are frequently referred to as cracks in fracture mechanics. The motion of a crack is 
determined by the stress field in the material under consideration. Based on a linear approximation , the classic equation of motion for a crack 
is beautifully described by the so-called linear elastic fracture mechanic s (LEFM) iFreund II 1998b. LEFM pre dicts that cracks are inertialess 
if the medium is unbounded, which has been experimentally verified by Goldman. Livne. & Fineberg 1 20ld) . 

On the other hand, if the medium is finite in size, Goldman et al. show that the motion of a crack is influenced by the boundaries of 
the medium so as to acquire inertia. This is possible because when a crack moves, the elastic waves surrounding the crack were reflected by 
the boundaries and have an influence on its motion. In other words, a moving crack interacts with its prior history. Goldman et al. further 
demonstrate that the inertia of a crack increases with crack velocity and becomes effectively infinite as the crack's velocity approaches a 
limiting velocity, cr, viz. the speed of Rayleigh wave in this medium. This increase of inertia parallels the mass divergence of particles in 
special relativity. As a result, the speed of the proposed S field should be equal to the speed of light. 

Then how a particle acquire inertia with the existence of external gravitational field? If we treat gravitational field as "boundaries" of 
the medium as for the case of crack, then a particle will acquire inertia in a similar way. The S field surrounding the moving particle will 
be reflected by gravitational field so as to influence the motion of the particle. If gravitational field vanishes, the particle will be inertialess 
because the particle is moving in a unbounded "medium". The stronger the gravitational field, the more the S field is reflected so that the 
particle's inertia is larger. When the gravitational field is strong enough that all S field is reflected, the particle's inertia will not increase and 
stays constant. 

What happens if the gravitational field vanishes and particles are therefore inertialess? Do they move at a velocity which is larger than 
the speed of light? As we said before, the S field should be at a velocity that is equal to the speed of light. When the inertialess particles' 
velocity approaches the speed of light, they will be strongly interacted by the S field they just transmitted so that they acquire inertia. In this 
limit the particles all behave like photons. Therefore we conclude that the speed of light is still the limiting speed. 

It should be stressed that inertia is only indirectly dependent on gravitational field. What gives a particle inertia is the interaction of the 
particle with the S field it transmitted previously. This point should be borne in mind when we derive gravitational field equation in Section 
231 

By these three assumptions, we shall derive the MOND relation ([2]) in next section. 



2.2 Results 

Now we consider one particle's motion under the gravitational interaction of a massive object at large separation. The particle moves inward 
by converting potential energy to kinetic energy: F = dp/dt = (dmi/dt) v + mi (dv/dt), where p = miv and mi, of course, is its 
inertial mass, as usual. Because F = — GMm g r /r 3 , where M is the gravitational mass of the massive object, we get 

dmi GMm K 

mm H — v = 7, — r. (3) 

at r 6 

For radial motion, we have 

, , 2 GMm R 
mivv + niiv — — -, (4) 

here the primes denote derivatives relative to the distance r. Substituting u for 1/r, the above equation is reduced to 

dv dmi 2 

miv 1 v — GMm g . (5) 

du du 

To move further on, we have to figure out, under the guide of the third assumption proposed in above section, how mi is shaped by external 
field. In general, based on the assumption that mi oc m g , an inertial mass function of the form 

mi = v m s (6) 

is expected, where v (x) is a function of jjn /oi only. As is well known, mi is a constant in the strong field case, and by the third assumption, 
mi — > in the weak-field extreme. We are interested in the weak-field extreme case. A plausible assumption for the weak field case is 
mi oc (<7N) a m g , or 
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mi (g N < ai) = J mg ' ( - 7 - ) 

where (/n = GM/r 2 — GMu' 2 , the Newtonian gravitational field. Here we introduce a new constant, Oi, which is related to Milgrom's 
constant do, as can be seen in the next section. In order to smoothly bridge the strong field case and weak field case, we expect < a < 1, 
which means mi increases rapidly when i?n <C ai but ceases to climb up when i?n — ax. The index a has to be fixed phenomenologically. 
As can be easily checked, if we set a = 1/2, the desired MOND behaviour is recovered, i.e. 

/ \ 1/2 

mi (g N < oi) = ( — j m g , (8) 



where 



ax 



oi = vo/GM, (9) 
and wo the particle's asymptotic velocity. Because v (x) « 1 when x 3> 1, a simple but quite plausible assumption is 

for all x. It should be stressed that, in this model, mi/m g is influenced only by external fields, not by any other factors, including the particle's 
dynamical quantities, e.g. its velocity. That is to say, mi is a true scalar. Because of this behaviour of inertial mass, the particle's Lagrangian 
can be expressed as equation (|13|) . 

If the particle is on circular orbits, therefore a constant inertial mass, the equation of motion, by yj, is reduced to 

m gflN = mxa, 

which reads 

a = (g N ai) 1/2 , 

i.e. the recovery of equation (|2|) if a\ = ao is recognised. This equation, however, cannot be applied to any other orbits other than circular 
ones due to the variation of the particle's inertia. In general, equation ([3| should be used. This indicates that the original MOND prescription 
is most suitable to describe the dynamics of spiral galaxies where internal motion is almost perfectly circular. 

Now let us consider how to escape the gravity of a massive object. Substituting equation JHJ into equation J5j we find, for a radial 
motion, 

UV- hi) =V . (11) 

du 

This equation tells us that if the particle is moving away from the massive object and has a velocity v = vo, then dv /du — and the particle 
keeps moving at a constant velocity. Therefore the particle is able to escape the gravitational pull of a massive object by itself. 

If v < vq, then dv/du > and the particle decelerates until reach a maximal distance. Of particular interest is the case of v > vq, 
where dv/du < and the particle keeps accelerating and the acceleration is increasing. It is straightforward to check that its velocity satisfies 
the following equation 

t 2 2\l/2 . 
[v — Vq) = Ar. 

Obviously, if no other field's disturbance, the velocity will increase infinitely so that v vo and 

v = Ar. (12) 

This is the direct result of the decreasing inertia while the particle moves away from the gravitational field. This equation is formally same as 
Hubble's law. But please bear in mind that Hubble's law describes velocity field, equation (|12[) . on the other hand, describes one particular 
particle's velocity change with distance, we shall defer the discussion of this equation to Section [2~6l 

It is clear that, since the particle's mass is dependent on the strength of gravity, gravity's apparent effect is not always attractive, but 
takes on different aspects dependent on the particle's motion state. If the particle is bound to the massive object, the gravity's apparent effect 
is attraction. On the other hand, if the particle is in an unbound state, the gravity's effect is repulsion, as indicated by (|12|l . the particle's 
velocity is not decreasing but increasing as it moves away from the gravitational field. 
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2.3 Lagrangian Formalism 

In above section we just assume that inertial mass depends on the external field and then extend Newtonian dynamics in a minimum manner. 
To formulate a self-consistent theory we need to develop a set of Lagrangian formulae. As usual, we write down the Lagrangian: 

1 2 

L = —miv — m g (p, (13) 

where <f> is the scalar gravitational field that is determined by Poisson equation V 2 <^> = AirGp with p the gravitational mass densit)Q. The 
Euler-Lagrange equations of motion 



dL 
At \dxi 



dL 

dxi 



= 



^ = -m g 'V<t>+ ^v 2 Vmi. (14) 

Compared with equation ([3]) we find that equation (|14[) contains an additional term iw 2 Vmi, stemming from the fact that inertial mass 
depends on external field, so is a function of position. It turns out that this term plays a key rule in accounting for the dynamics of dwarfs, 
spiral galaxies and clusters in a consistent and self-contained way. Applying p = miv to the left hand side of equation (|14[) gives: 

1 2 

mia = m g gN + — v Vrrer — (Vmi • v) v. (15) 
For the deep-MOND case, substituting equation ([8]) into above equation yields 

— gN + 7 — V# N -- v. (16) 

This equation is the general motion equation for the deep-MOND regime, including the non-spherically symmetric systems. It is apparent 
from equation (1151 ) that in a modified inertia theory, unlike modified gravity, the acceleration experienced by particle depends not only on 
position, but on the velocity as wel l. This seems a drawback of this kind of theory, but on the other hand this could be an unparalleled 
advantage (see, e.g., iMilgroimlbooi) . 



As above, let us first consider the circular motion. In this case the above equation, after applying equation ([8]), reduces to: 



a =^ )1/2 { 1+1 27m=)- (17) 



From this equation we can easily find the stable circular velocity 



v c = (4GMai) 1/4 . (18) 

However, only the sufficiently virialised systems can attain to this high circular velocity. For those systems less virialised, the objects on the 
outskirts of the systems will follow quasi-circular orbits and gradually wind up during the inward migration. This indicates that Milgrom's 
constant is actually slowly varying according to 



do 



a x 11 + 



(19) 



that is to say, ao is varying in a narrow range at ^ do 4ai and older systems tend to have higher values of ao. 

To find the value of ai, an easy but reliable way is to carefully select a large enough sa mple of w ell studied galaxies so that t he scatter 
of ap is m oderate and MOND works qu ite well for this sample, we select the galaxies from lBegeman. Broeils & Sanders I Jl991 ), Sanders 



and lSanders & Verheiien 1 1 119981) . The reasons for selecting these galaxies are three-fold. First, these galaxies are among the highest 
quality of observational data and it was demonstrated that MOND can account for the data with a relatively high precision. Secondly, these 
galaxies are all spiral galaxies for which MOND is most successful. In addition, as indicated above, ao evolves with galaxies. To reduce the 
scatter of ao, we should select galaxies with comparable properties. Finally, these three samples contain several galaxies, e.g. NGC 2841, that 
are quite controversial within the context of M OND. It is quite desirable if a varying ao can settle down these issues. The resulting sample 
contains 63 spiral galaxies, but 9 galaxies from I Sanders & Verheiien I j 19981) are eliminated from the list adopted to calculate ao because of 
the reason presented below. 



1 As stated above, when we vary 0, we should get Possion equation because mj is not directly dependent on gravitational field. 
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In Table 1 the calculated values of ao are listed for every galaxy in the sample, ax is found by requiring the average value of ao = 
1.21 x l(T 8 cms- 2 : 

ai = 0.667 x 10~ 8 cms" 2 . (20) 

Several points should be mentioned. The accelerations at the last measured points of the rotation curves (ai mp ) of NGC 3949, NGC 3953, 
NGC 4085, UGC 6973 are in the Newtonian regime and therefore eliminated from the calculation of oi. In light of the found value of ax, we 
expect that the galaxies with ai mp comparable to or larger than ax will be quite dynamically different from other galaxies, these galaxies are 
NGC 3877, NGC 3972, NGC 4051, NGC 4217, NGC 4389, with a imp > 0.7 x 10" 8 cms" 2 . As a result, these five galaxies are eliminated 
from the calculation of ax. But it should be mentioned that, if these five galaxies are included in the galaxies responsible for calculating Oi, 
ax would be slightly smaller. By this selection criterion, the resulting scatter of ao, ±0.30, around the mean value is moderate, as desired. 



2.4 Comments on Selected Individual Galaxies 



A quick glance of the result for ao suggests that large and high surface brightness galaxies tend to have a high ao. Several comments are 

presented below for those values significantly deviated from the standard value: 

NGC 2841. - This is the most controversial spiral galaxy for MOND. iBegeman. Broeils & Sanders" 1 1991 ) found that, if a Hubble 
distance (9.5 Mpc) was used, MOND obviously fails to account for the rotation curve. On the other hand, if the Tully-Fisher distance 
was used, MOND is in well agree ment with the observ ed data. However, the Tully-Fisher distance is twice as the Hubble distance, a seemly 
unacceptable result. Subsequentlv. lBottema et al. I (12002!) found a distance of 14.1 Mpc to NGC 2841 based on Cepheid method. This distance 
goes half way between the Hubble distance and Tully-Fisher distance. As a result, the situation is alleviated somewhat, but the predicted curve 
still deviates systematically from the observed data. Table 1 gives this galaxy's ao a value of 1.48 x 10 -8 cm s -2 , the largest in this sample if 
not include the values not adopted for calculating at, which is more than 22% larger than the standard value. This value will further alleviate 
the situation and make MOND in principle compatible with the data. But it should be pointed out that even with this high value, the data still 
can not be comfortably matched by the prediction. Here is a caveat. We find that, for this galaxy, ai mp = 0.66 x 10 -8 cm s -2 ~ ai. Because 
ax is the transition acceleration beyond which the modified inertia transits to Newtonian inertia, we expect the value ao calculated in this way 
is only qualitatively correct. If otherwise ai mp <C ai, the value listed in Table 1 should be exact. But NGC 2841 evolved beyond this stage 
because it has a much higher ai mp . Quantitatively, a value of 1.87 x 10 -8 cms -2 for ao will do the work if the Cepheid distance is used. If 



we adopt equation (|10|) as the inertial mass dependence on external gravitational field, we find a value of 1.86 x 10 cms - (~ 2.77ai) 
for ao, which is just what is needed for this galaxy to bring MOND prediction in accordance with the observation. As in the general case, 
when ai mp ~ oi, the actual ao depends on the ratio of rotation velocity to its stable circular velocity. When rotation velocity approaches its 



stable circular velocity, ao will have a value of 2.45 x 10 cms (~ 3.66ai). This implies 2.77ai ^ ao ^ 3.66ai when ai„ 



a x . We 



find that the maximum value (3.66ai) for ao when ai mp is in the vicinity of ai is less than the corresponding value (4ai) when ai mp <C ax, 
indicating that the Newtonian dynamics begin to take over when acceleration enters the Newtonian regime. The oscillating rotation curve, 
rather than a perfect flat rotation curve, of this galaxy on the outskirts may have some implication of this subtlety. 

DDO 154. - A value of 0.9 x 10 -8 cms -2 , significantly less than the average value, is given for this dwarf and gas-rich galaxy. 
Nonetheles s, this low value can be larg ely attributed to the low rotation velocity at the last measured points. This is a q uite controv ersial case 



for MOND. Milgrom & Braun 1988D demonstrated that this galaxy is among the most acute t est of MOND. Sh ortly, Lake i 19891) , however, 



showed that a substantially small value of ao is necessary to account for the observed curve. Milgrom] jl99lh disputed Lake's conclusion, 
arguing that a different distance to this galaxy be more reliable. But if the ideas presented in this paper are correct, a small value of ao is 
inevitable because of the declining nature of the rotation curve at the outermost points. The internal acceleration of DDO 154 is everywhere 
less than ao so that MOND should play a dominant role. But MOND cannot account for the declining nature of the rotation curve. In the 
spirit of this paper, however, the declining nature of gas-rich dwarf galaxies can be naturally accounted for. The main point is that gas-rich 
galaxy is newly formed object with a less relaxed outer part and a fairly relaxed inner part. Consequently, the inner part should have a higher 
value of ao than its outer part. This indicates that ao is varying not only from galaxy to galaxy, but also within galaxies. For sufficiently 
relaxed galaxies, ao could be treated as a constant without loss of accuracy, but for gas-rich galaxies, the variation of ao is maximised. 

NGC 3893. - Table 1 gives ao a value of 1.47 x 10 -8 cms -2 . This high value is largely the result of the high ai mp . But because this 
galaxy has a disturbed velocity field, the result should not be taken seriously. 



NGC 3992. - Table 1 gives ao a value of 1.48 x 10 cms - or 1.86 x 10 cms - if equation (|10l) is used, an apparently not 
necessary adju stment for the standa r d valu e if the good agreement between predicted and observed rotation curve is noticed. However, as 
mentioned by Sanders & Verheiien 1 1998), the fitted value of mass-to-light ratio of the stellar disk is unusually large compared with other 
values in that sample. If, however, an ao with above value is used, the near-infrared mass-to-light ratio will be reduced by a factor of 1.22 or 
1.54, a value compatible to other galaxies in the sample. 

NGC 4010. ~ Table 1 gives ao a value of 1.44 x 10 -8 cms -2 and even higher if equation (|10l) is used. This high value of ao seems 
unexpected for such a low surface brightness gal axy. However, this galaxy has a very extended disk and the scale length is much larger 
compared with other galaxies in its mother cluster ITullv & Verheiien II 19970 . As a result, this galaxy has a much larger internal acceleration, 
so is ao- This galaxy is an example that ao does not actually directly depend on surface brightness. Since young galaxies with outstanding 
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star formation activities tend to have a high surface brightness even if its internal acceleration is low. On the other hand, some old galaxies 
devoid of star formation tend to have a low surface brightness because of their old star population. In addition, we find that the rotation curve 
of this galaxy begins to decline at large radii, anothe r example that ao is varying within one galaxy. 

Recentlv lSwaters. Sanders & McGaugh 1 020101) study a sample of 27 dwarf and low surface brightness galaxies within the framework of 
MOND and find that there are some systematics that cannot be reconciled easily: MOND curve predicts higher rotation velocities in the outer 
regions and lower in the central regions for low surface brightness dwarfs; higher rotation velocities in the central regions and lower in the 
outer regions for high surface brightness galaxies (UGC 5721, UGC 7323, UGC 7399, UGC 7603, UGC 8490). This intriguing phenomenon 
forces Swaters et al. to try a fit with ao as a free parameter. By allowing ao vary, Swaters et al. find a trend that lower surface brightness 
galaxies tend to have lower ao. This is just what equation U9I) expected if lower surface brightness galaxies tend to be less relaxed. If equation 
( 1191 ) is app lied to this sample of systems, an aver age value of 1.06 x 10 -8 cms -2 is found, which is significantly lower than the commonly 
used value. Swaters. Sanders & McGaugh 1 20ld) . however, find an even lower average value of 0.7 x 10~ 8 cms -2 for ao in this sample, 
a value that is just above oi. Though the calculated value of ao is significantly lower than the commonly used value, it is still significantly 
higher than the fitted average value in this sample. Furthermore, the fitted values for the galaxies in this sample have an enormously large 
scatter that is completely out of the reach for (1191) to account for. This enormous scatter, on the one hand, is a consequence of the internal 
dynamics of the systems; on the other hand, is a consequence of observational uncertainties. If such a large scatter of ao is necessary for this 
sample, then it will be a challenge for the ideas presented in this paper. 



2.5 Dynamical Mass of Clusters of Galaxies 



It has been quite controversial with the applicability of MOND to clusters. iThe & Whi te ( 1988) first noted that MOND is inconsistent with 



the dyn amics of Coma cluster unless the MOND acceleration parameter ao is 4 times larger than the value inferred from spiral galaxy rotation 



curves. SandersJ i ll 9941) . based on a small sample of X-ray-emitting clusters, concluded that MOND mass is consistent with the detected gas 



mass. Subsequently, based on a much larger sample o f X-ray-emitting clusters, Sanders found that MOND over-predicted the mass by a 
factor of 2 than the detected gas mass I Sanders] [l999l) . In light of this finding, the mass over-prediction was already implied bv lSanders 
1 1994) though with a less statistical significance. 



This m ass over-predicti on flaws MOND since MOND was particularly speculated to eliminate non-baryonic dark matter. To rem- 
edy MOND, I Sanders I d2003l) suggested that neutrinos, with mass of 2 eV, is responsible for the mass discrepancy in rich clusters. This 
mass of neutrinos, having a neg ligible effect on the galaxy scale, could have a significant dynamical effect on cluster scale. Unfortunately, 
Angus. Famaev & Buote I f2008T) showed that neutrinos are incompatible with the observed mass distribution within X-ray bright groups and 
clusters. C onsequent l y, Ang us et al. proposed that sterile neutrinos would close the cluster mass problem. This conclusion contradicted the 
finding by Milg roml jl998h who fou nd that MOND is compatible with galaxy groups. Unfortunately, by studying the gravitational lensing 

of clusters. iNataraj an & Zhao I OOP 81) conclude that even sterile neutrinos are far from enough to close the mass budget of cluster s. 

19991) . the 



However, if the inertia is allowed to vary, as presented above, no dark matter is needed. Of the 93 clusters studied bv lSanders 



internal accelerations are of ~ 0.6 x 10 cm s~ , that is, in the vicinity of at- Section l2.4l shows that for systems with internal accelerations 
comparable to at, a much higher effective ao should be used. It can be checked that for these 93 clusters, ao approaches the maximum value 
3.66ai (2.02 x 1.21 x 10~ 8 cm s -2 ) . This high value is just enough to eliminate the mass discrepancy in bright clusters. In this respect, we 
say that the dynamics of bright clusters are quite similar to that of NGC 2841. This also implies that not all clusters would have the mass 
discrepancy problem within the framework of MOND. Those clusters with much lower internal accelerations will have similar dynamics as 
spiral galaxies. 

But the mass discrepancy problem is not completely closed because in some cases MOND over-predicts masses by a factor of 3 or even 
higher. If this phenomenon is ubiquitous and no conventional interpretations, e.g. dark baryons or neutrinos, are reliable, then there really 
exists a dark matter problem for MOND. 



2.6 Unbound Systems 

After a discussion of circular motion, let us now turn to radial motion. In this case equation (|15l) reduces to 

a = (aigr-j) 1 ^ 2 ( 1 — — . V — I . 

We find the conclusion is not changed much for radial motion by an application of Lagrangian formula. The only difference is a change of 
the asymptotic velocity, defined by equation l|18|) other than the value vo determined by equation (|9|. 

Equation (|12[) reminds us of the Hubble's law. But does it really imply Hubble's law. By its low magnitude (Ho = 72 km s -1 Mpc -1 ) 
compared to the typical galaxies' asymptotic rotation velocities (~ 200 km s -1 ) and their typical sizes (< 100 kpc), if it is really Hubble's 
law, we expect that Hubble's law is not the result of mass inhomogeneity on the galactic scale. On contrary, Hubble flow is a result of mass 
inhomogeneity on cosmic scale. Obviously, if the mass distribution is completely homogeneous, galaxies would not feel the existence of 
an external gravitational field, resulting in zero inertial masses and therefore random motions. But actually every galaxy is embedded in an 
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ubiquitous gravitational field. A fluctuation of mass density on cosmic scale gives rise to the small gravity on large scale, those cosmic objects 
unbound by this inhomogeneity of mass distribution acceleratingly expand. If 8M is the mass inhomogeneity caused by a density perturbation 
A = Sp/p on the cosmic scale, we expect the Hubble recession velocity exceeds 5vo = (AGSMaj) 1 ^ 4 so that an accelerating expansion 
is guaranteed. This indicates that the accelerating expansion of the Universe can be most readily observed on scales that the clustering of 
galaxies becomes negligible. The two-point correlation function of galaxies shows that these scales are r ~ 5/i _1 Mpc. So the accelerating 
expansion of the Universe should be free of peculiar contamination beyond these scales. The accelerating expansion of the Universe seems to 
imply an increase of the universe energy, or in modern parlance, the existence of dark energy. But above analysis shows that the accelerating 
expansion of Universe does not violate conservation of energy because equation ()12|) is derived based on the conservation of energy. All 
difference is the dependence of inertia on gravity. Although galaxies' velocities are becoming larger and larger, their inertial masses are 
dropping dramatically, resulting in the conservation of total energy. 

Although we recover MOND's relation ([2]) in the case of circular orbits, the ideas presented above are quite different from that of 
MOND. In MOND, gravity's effect is always attractive, the moving particles in the gravitational field are all feeling the same acceleration as 
long as they locate at the same place in the field. In fact MOND, including the relativistic MOND theory proposed by Bekenstein, complies 
with weak equivalence principle. This is however not the case for the ideas presented here. Just as we already find, gravity is attractive 
in the case of bound state, and repulsive in an unbound configuration. Even for particles in the same type of state, say bound state, their 
accelerations could be quite different because of their different radial velocities, which determine at what rate the moving particles' inertial 
masses change. As a consequence, all particles moving from large distance towards the massive object tend to converge to the same velocity, 
i.e. the asymptotic velocity vo = (4GA/ai) 1 ^ 4 . If the particles' velocities are less than vo, the attraction effect of the gravity increases their 
velocities. On the other hand, if the particles' velocities are larger than vo, the repulsion effect of the gravity decrease their velocities. Only 
in the Newtonian regime is the gravity always attractive. 

In addition, since the gravitational potential is not infinite, opposed to the original MOND, the moving body can escape the gravitational 
pull easily without resorting to other nearby massive galaxies' assist as described bv lFamaev. Bruneton & Zhao I d2007l) . This is yet another 
advantage of the changing inertia proposal over the traditional MOND. In many cases, by finely tuning the parameters, the predictions 
of CDM model are barely discernible from that of MOND. But globular clusters are outliers, which are believed contain no dark matter. 
Consequently iBaumgardt. Grebel & Kroupal ( 120051) suggested a test of gravitational theories by studying velocity disper sion of glo bular 



clusters. Recently these studies are carried o ut by several groups, finding an inconsistency between MOND and observations (Moffat & Toth 



2008l : lHaghi et al. l2009l : ISollima & Nipoti hoifj) . The problem is that MOND systematically overpredicts the velocity dispersions. This can 



be understood as follows: Because of the infinite escape velocity intrinsic in MOND, at any given radius the velocity distribution should be 
large. This is no longer the case for the changing inertia proposal since at any radius the escape velocity is finite. 



3 DISCUSSION 



3.1 Indications for structure formation and evolution 



The cosm ology and s tructure formation in the context of MOND were first studied bv lFeltenl a 1984i) and then extended by Sanders ( 1998), 
Sanders I ( 120011) and lStachniewicz & Kutsc hera (2001). T he main findings of the se studies are that MOND would boost the structure for- 
mation even with a very low matter density. Ivan den Bosch & Dalca nton (2000) semianalytically studied the formation of disk galaxies 



and found that dark matter and MOND models are indistinguishable. IStachniewicz & Kutschera (2005) found that the dark ag es would 



end at z ~ 30, in good agreement with the WMAP results. More recent studies (Sanders 



20081 ; iMalekiani. Rahvar & Haghi H2OO9I) of 



structure formation generally str engthen these early conclusions. When confronting MOND with the formation of galaxy clusters, however, 
Nusser & Pointecouteau (2006) found that the simulated MOND clusters are significantly denser than the observed ones. It means that 



MOND, as a gravitational field, is so strong that it over-evolves the structure on cluster scale. 

If the ideas presented in this paper are correct, we find that the changing inertia postulate boosts the structure formation but the evolution 
is only moderately boosted. This is because when the particles begin to fall into the potential well of the density perturbation, thanks to their 
low initial velocities, therefore small variation of their inertial masses in unit time duo to small variation of the potential well in which they 
are moving, the accelerations of the particles are almost the same as that predicted by the standard MOND. This indicates that the timescale of 
initial formation is barely different from that of standard MOND. But when the particles gain enough velocities, owing to the rapid variation 
of their inertial masses, their accelerations are quite small compared with the ones predicted by the standard MOND. There fore the structure 



dynam ical evolution is deferred. This effect will make the structure less dense than that predicted by MOND. Ivan den Bosch & Dalcanton 



(2000) also found that stellar feedback is needed to reproduce the lack of high surface brightness dwarf galaxies within the context of MOND. 
However, the changing inertia proposal will defer the evolution of structure so that central surface brightness could be significantly lower. 

This line of arguments also indicates that older objects in a system should be kinematically hotter than their younger counterparts. I t 
is well known that there exists an age-velocity relation (AVR) in the disk of spiral galaxies fcf. iMeusinger. Stecklum & Reimannlll99lh . 
The ideas presented in this paper indicate that there should be a rapid rise of the velocity disp ersion of disk stars and the n saturate when 



they a re old enough. Unfortunately, the kinematical evolution of the disk is still controversial. Edvar dsson et al. ( 1993) and Soubiran et al. . 
( 2008 ) found that there is a rapid rise of velocity dispersion for the first 2 Gyr and then saturate for old er stars. iNordstrom et al. 1 12004|) and 



Rocha-Pinto et al. 1 120041) . on the other hand, found a steady increase in the velocity dispersion with age. lHolmberg. Nordstrom & Andersen 
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(2007), with an improvement on T e g and [Fe/H] calibrations for early F stars, found a steady increase in velocity dispersion with age and a 
tentative saturate at age 10 Gyr. 

The age-velocity relation is usually attributed to disk heating. Several heating mechanisms, e.g. scattering by giant molecular clouds, 
spiral structures, and accretion, have been pr oposed to account for th is relation. But it is growingly evident that these mechanisms are 
insufficient to account for the observations (cf. Rocha-Pinto et al.1l2004 ). The ideas presented in this paper may shed new light on this issue 
because they will boost the velocity dispersion. 



3.2 More words on inertia 

When we argue for a dynamical origin of inertia, we mentioned cracks and the acquisition of their inertia as an interaction with their prior 
histories by introducing boundaries for the medium. Here we present evidence that why cracks are of particular interest. Cracks, the line 
defects in material, have their siblings in crystal, viz. dislocations and disclinations. Since the development of defect gauge field theory 
in condensed matter physics, it has been more and more clear that Ei nstein's gravitat ional theory is a special case of the gauge theory for 
defected crystal in the sense that general relativity is free of torsion (cf. lKleinert 1990). The geometrical study of crystal finds that a crystal 
filled with dislocations and disclinations has the same geometric properties as an affine space with torsion and curvature, respectively. 

The Einstein curvature tenser, Gp V , in its three-dimensional version, is equivalent to the disclination density tenser, 0;j. As a result, 
the Einstein field e quation, G uv = — 87rGT M „, indicates that the baryons are just disclinations in the "world crystal" with Plank length as 
the lattice constant ( Kleinert 1990). Then why do we not observe a space with torsion? Torsion can only be a result because of the existence 
of dislocations. However, dislocation would make the space discontinuous and cannot exist macroscopically because it is unable to maintain 
stress equilibrium by itself. This is, however, not to say that dislocations do not exist microscopically. Actually the dislocation density, a^, 
is related to disclination density by diOHj — —Sjki&kh That is to say, several dislocations could end up as a disclination as seen by a distant 
observer. Therefore dislocations are akin to quarks in the way that dislocations and quarks are always confined within disclinations and 
baryons, respectively. If quarks are indeed dislocation defects in world crystal, then we should expect that space is not torsion free on the 
quark scale. Furthermore, the conservation laws for the angular momentum and energy-mom entum density in general relativity are just the 
respective conservation laws for dislocation and disclination density in defected world crystal i Kleinert 1 199fj|) . 



This indicates that the argument for the way particles acquire inertia in a similar way in which a crack acquires inertia in a bounded 
medium is not just a trick but a facet of reality. 



4 CONCLUSIONS 

MOND has been well established as a serious alternative to the standard dark matter model for about three decades owing to its remarkable 
success in accounting for the dynamics of a variety of galaxies with quite different luminosities, morphologies. However, when confronting 
with dwarfs and clusters, MOND is controversial. There exists cumulated evidence that the acceleration parameter ao, assumed to be univer- 
sal, is varying in the sense that low surface brightness galaxies tend to have low ao- For bright clusters, a factor of 2 mass over-prediction is 
well established. Again, this over-prediction can be accommodated by an ao 2 times larger than the usually adopted value. 

With above observational evidence, a motivation to theoretically account for the phenomena is prompt. This paper postulates that the 
MOND phenomenology can be accounted for by three assumptions: 1) Gravitational mass is conserved; 2) Inverse-square law is applicable 
at large distance; 3) Inertial mass depends on external fields. The first two assumptions are quite general and should not find objection. The 
third assumption, given by equation ® appropriate at the low acceleration regime, is key to reproduce the phenomenology of MOND. It is 
found that the inertia modified in this way exactly recover the formulae suitable to MOND in the circular motion case. 

By a consideration of Lagrangian formulae, however, it is found that the usual relation (|2j is replaced by (|17|) . A comparison of equations 
([2]) and (|17|) shows that ao, a proposed universal constant, is actually varying in a narrow range, ai ^ ao ^ 4ai, if the internal accelerations 
of the systems in question (or the external field) are much lower than ai. A scrutiny of equations ()15[) and (|10l) indicates that the effective ao 
not only depends on the orbital velocity but also on the external field. This varying ao is just enough to eliminate the over-prediction factor 
of 2 found in bright clusters and the lower value in low surface brightness galaxies. Because of this nature of ao, ao not only varies from 
system to system, but also varies within one galaxy. This phenomenon can readily account for the decline of rotation curves found for many 
low surface brightness galaxies. 

Above statement suits only for circular motion. For radial motion, situation is quite different. It turns out that there exists a critical 
velocity, v c defined by equation l|18|) . for every galaxy. If the object in radial motion is moving slower than this critical velocity, gravity 
appears to be attractive, as usual. But for a motion faster than this critical velocity, gravity appears to be repulsive in the sense that the 
object begins to accelerate. This immediately reminds us of the accelerating expansion of the universe. By a careful inspection of the cosmic 
parameters, if the modified inertia is correct, it is realised that the accelerating expansion of the universe must be the result of inhomogeneity 
of the universe on cosmic scales. 

In general, besides circular motion and radial motion, equations (|15|) and (|10|) should be used to account for the dynamics of systems 
in question. These two equations indicate that the traditional MOND prescription is only appropriate to describe circular motions. It is found 
that ao is not a fundamental constant, its value is an indication of the system's age and relaxation index. In general, high surface brightness 
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objects tend to have high value of do. But this is not appropriate for all systems, ao is actually determined by the orbital velocities and 
external gravitational fields. A more accurate value of ai should be determined by these two equations. But owing to the success of MOND 
in spiral systems, the true value of a\ should not deviate far from the value given by equation Q20[) . 
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Table 1. The sample se lect ed for the calculation of au C alculated values of qp are listed on column 2. Refe rence: (l) lBegeman. Broeils & Sanders 



1 199 lh : (2) ISanderslfl996h : (3) ISanders & Verheijen I ( fl998h . Nine galaxies from lSanders & Verheijen I Jl998h , NGC 3877, NGC 3949, NGC 
3953, NGC 3972, NGC 4051, NGC 4085, NGC 4217, NGC 4389, UGC 6973, are eliminated from the calculation of aj because of the rea- 
son presented in text. The data for UGC 6446 are taken from lSwaters. Sanders & McGaugh U2010I) to utilise the most recent update, but the 
differences are small. 
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Galaxy ao (10 8 cms 2 ) ref 

NGC 2403 1.20 1 

NGC 2841 1.48 1 

NGC 2903 1.32 1 

NGC 3109 1.01 1 

NGC 3198 1.11 1 

NGC 6503 1.06 1 

NGC 7331 1.29 1 

NGC 1560 1.12 1 

UGC 2259 1.24 1 

DDO 154 0.90 1 

DDO 170 0.99 1 



UGC 2885 1.28 2 

UGC 5533 1.17 2 

NGC 6674 1.16 2 

NGC 5907 1.34 2 

NGC 2998 1.20 2 

NGC 801 1.12 2 

NGC 5731 1.24 2 

NGC 5033 1.24 2 

NGC 3521 1.25 2 

NGC 2683 1.31 2 

NGC 6946 1.16 2 

UGC 128 1.00 2 

NGC 1003 0.98 2 

NGC 247 1.22 2 

M33 1.32 2 

NGC 7793 1.36 2 

NGC 300 1.09 2 

NGC 5585 1.10 2 

NGC 2915 1.06 2 

NGC 55 1.15 2 

IC2574 1.05 2 

DDO 168 1.14 2 

NGC 3726 1.20 3 

NGC 3769 1.02 3 

NGC 3893 1.47 3 

NGC 3917 1.33 3 

NGC 3992 1.48 3 

NGC 4010 1.44 3 

NGC 4013 1.26 3 

NGC 4088 1.38 3 

NGC 4100 1.31 3 

NGC 4138 1.31 3 

NGC 4157 1.30 3 

NGC 4183 1.11 3 

UGC 6399 1.24 3 

UGC 6446 1.05 3 

UGC 6667 1.23 3 

UGC 6818 1.17 3 

UGC 6917 1.31 3 

UGC 6923 1.44 3 

UGC 6930 1.23 3 

UGC 6983 1.15 3 

UGC 7089 1.14 3 



